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ABSTRACT 

It is demonstrated that decimation of the one dimensional Ising model, with periodic 
boundary conditions, results in a non-linear renormalisation transformation for the 
couplings which can lead to chaotic behaviour when the couplings are complex. The 
recursion relation for the couplings under decimation is equivalent to the logistic map, 
or more generally the Mandelbrot map. In particular an imaginary external magnetic 
field can give chaotic trajectories in the space of couplings. The magnitude of the field 
must be greater than a minimum value which tends to zero as the critical point T = 
is approached, leading to a gap equation and associated critical exponent which are 
identical to those of the Lee- Yang edge singularity in one dimension. 



The renormalisation group has been developed into an immensely powerful tool for the 
analysis of physical theories near critical points and also for continuum field theories. 
There are by now various forms of "renormalisation group equation" which govern 
how physical amplitudes and couplings change under change of scale. Perhaps one 
of the most intuitively appealing is the version due to Wilson [1] , motivated by a 
suggestion of Kadanov [2] , involving "decimation" . 

In principle one can derive recursive formulae for the couplings of a theory, which 
dictate how they should change when the underlying lattice is decimated, so that 
the Hamiltonian involving the new couplings on the new lattice is the same as the 
Hamiltonian involving the old couplings on the old lattice, i.e. the partition function 
does not change under the simultaneous operations of decimation and redefinition of 
couplings. The recursive formulae for the couplings are in general non- linear (indeed 
they are not invertible, so the transformation involved here is not a group but a 
semi-group) . 

Non-linear recursive formulae are one of the central themes of study for chaos theory 
and one can pose the question - can the renormalisation transformation lead to chaotic 
behaviour in the space of couplings? This possibility has been investigated before, 
and answered in the affirmative using numerical calculations in some specific models 
for which exact recursion relations can be obtained, [3] [4] [5] [6] . Here a simple model 
(the one dimensional Ising model) will be analysed analytically and it will be shown 
that this model also exhibits chaotic behaviour in a surprisingly elegant manner. The 
analysis shows that the onset of chaotic behaviour appears to be associated with the 
second order phase transition at T = 0. It remains an open question as to whether 
this is a peculiarity of this model or is a more general feature. 

One severe problem in extracting general features is the paucity of models for which 
the recursion relations can be obtained exactly, and the one dimensional Ising model 
- because of its simplicity - is one example for which progress can be made. Never- 
theless, despite its simplicity, the results are startling enough to merit description. 

It will be shown that the onset of chaotic behaviour is brought about by extending the 
couplings of the theory to the complex plain. This is not a new idea in the analysis 
of such theories. Dyson [7] pointed out that one could learn something about the 
structure of Quantum Electrodynamics by considering imaginary electric charges, so 
that a = e'^/hc < 0. Such a theory must be intrinsically unstable, and so amplitudes 
cannot be analytic at a = 0, hence perturbation theory must diverge and expansions 
in a are, at best, asymptotic. These ideas have been further developed by making 
a complex and there is by now a whole literature on complex analyticity and Borel 
summability (e.g. [8] ). In statistical mechanics, extending the couplings to the 
complex plane is a key step in solving many two-dimensional models [9] and has led 
to some beautiful results concerning the analyticity of the partition function [10] . 



2 



In this paper yet another example of the fascination of complex variables will be 
exhibited - by allowing the couplings of the one dimensional Ising model to be com- 
plex, the recursive rcnormalisation transformations can become chaotic. To exhibit 
this phenomenon, some well known features of the one dimensional Ising model will 
be summarised and the recursion relations derived. It will then be shown that the 
recursion relation is nothing other than the logistic map, and chaos ensues! 

Consider the one dimensional Ising model on a periodic lattice of N sites [9]. The 
partition function is 



Zn = ^ exp 



N N 



(1) 



where K = and h — with J the spin coupling and H the external magnetic field, 
(periodic boundary conditions require crjv+i = cri). ZN{K,h) can be conveniently 
expressed in terms of the transfer matrix 



V 



V++ V+. 



(2) 



as Zn = TrV 



N 



Diagonalising V gives the eigenvalues 

A± = jcosh/i ± Vsinh^/i + e-^^} 

Thus 



N 



Zn — 



1+ 



N 



(3) 



(4) 



The recursive rcnormalisation transformation is well known for the 1-D Ising model 
[11] . It is obtained by asking: can one find new couplings K' and h' such that 

Zn (K', h') = A^Zn {K, h) (5) 

gives the same physical amplitudes? [A is a normalisation factor.) 

Equation (5) is easily satisfied by demanding 
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giving the recursive formulae: 



2h' _ 2h cosh(2K + h) 
cosh{2K-h) 



^j^i _ cosh(4K) + cosh(2/j,) 
2cosh^(/i) 

The normahsation factor A is unimportant for the present analysis. 



The combination sinh^(/i') = e^-'^ sinh^(/i) is a renormalisation transformation 
invariant. This means that the magnetisation per site, in the thermodynamic limit, 

dlnZN e^^smh{h) 
is invariant under the renormalisation transformation. 

All these facts about the one dimensional Ising model are well known [9] and are 
included only for completeness. 

It will now be shown that the recursion relations (7) and (8) are equivalent to the 
logistic map and, for certain (complex) values of couplings gives rise to chaotic be- 
haviour. 

Define 

m= l + e^^sinh2(/i) (10) 

which is a renormalisation transformation invariant, m — m' . It is now only necessary 
to consider one of equations (8) and (7) as the existence of the invariant, m, makes 
one of them redundant. 

Eliminating h from (8) using (10) gives 

4[(e4^-l)+m]" ^ ' 

Now replace K with a new variable 

^ = -(?^ (12) 
with -oo < X < Q for m > and K > Q. The recursion relation (11) now becomes 

x' = Ax{\ - x) (13) 
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which is the logistic map. 



For < a; < 1, the recursion relation (13) leads to chaotic behaviour, as is easily seen 
by defining x = sin {nil;), < V < I (see e.g. [12] ) giving 

sin(7rV'') = sin(27rV') (14) 

Writing t(j in binary form, we see that the iterative map merely shifts all bits one 
step to the left and throws away the integral part, leaving the fractional part behind. 
For an initial value of which is rational this will lead to a periodic orbit, but for a 
starting value of ip which is irrational, the process never repeats and jumps around 
chaotically. Since the irrationals have a greater cardinality than the rationals, almost 
all initial values lead to chaotic motion. 

For real values of the couplings, m > 1 and -oo < a; < 0. Chaotic trajectories require 

m= l + e^-^sinh2(/i) < 0. (15) 

For example, if > is real and h is pure imaginary {h = i9), then m < for 
sin^ e > e-^^ and x = ^^"IIk'i"^ lies between and 1 for e^^ > 

The region of chaotic flow is shown in Figure 1. Note that not all points above the 
line X — lead to chaotic trajectories, only those with irrational '0. Indeed there is 
an infinite number of periodic trajectories as well as chaotic ones. For example there 
are lines of unstable fixed points (period one) at x = and x = 3/4 and the values 
X = (5 ± vTo) / 8 give orbits of period two , etc. 

For finite K there is a gap and a small imaginary magnetic field is not sufficient to 
induce chaos, but as K increases (T 0) this gap reduces to zero. Following Baxter 
[9], define t — e~^^ with t — > being the critical point, then near t = the line 
separating chaotic from regular flow is given by 

t'^r^e^ Qr^t. (16) 

If we define a critical exponent. A, such that the critical value of is 

e^t^, (17) 

then one obtains A = 1 for the one dimensional Ising model, which is just the critical 
exponent for the Lee- Yang edge singularity in this model. It is not difficult to see 
that the critical line e^^ siv? ^ = 1 is related to the Lee- Yang zeros of the partition 
function in the complex /j,-plane. The equation Zn{K, h) = has N roots, for real 
K, 

Zat = (A+)^ + (A_)^ = ^ A+ = e'^A_, (18) 
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where —N < q < N is odd. Using the exphcit form of the eigenvalues, (3), this leads 
to 

cos(l^) ^J e-^^ + smh\h) = ^sin(|^) cosh(/.). (19) 
This equation can be re-arranged to give 

cos(^,) = A/r~^cos (1^) , (20) 

where t = e~^^ as before and 9q are the roots in the rotated complex /i-plane, 
h — 19. Since < t < 1 the N values of 9q are all real and they lie in the range 
t < I sin(6'q)| < 1, which is precisely the region above the critical line in figure 1. In 
the thermodynamic limit {N — > oo) the lowest zeros are at sm.{9q) = t which is the 
critical line. Thus the critical line co-incides with the rightmost Lee- Yang zeros in 
the complex activity plane in the thermodynamic limit. 

One can obtain further insights by allowing K to become complex. Define x = — f + 1 
in equation (13) to give 

z' = z'^ - 2. (21) 

This is the Mandelbrot map z' = z'^ + c for complex z, with c = —2. One can have 
divergence or convergence depending on c and the initial choice of z. The values of 
c for which the iterates of the starting point z = stay bounded is the Mandelbrot 
set and clearly c = —2 is an element of this set. The Julia set for a given value of c 
is the set of points in the complex 2;-plane which stay within a bounded region upon 
repeated iteration of the Mandelbrot map (strictly speaking this is the filled in Julia 
set Jc - the Julia set J is actually the boundary of this set). The set J is generated by 
the inverse set of the unstable fixed points. For c = — 2 these are z — —1 and z = 2 
and so the earlier analysis of the logistic map tells us that the inverse iterates generate 
a dense set of points in the segment of the real axis lying between —2 and +2 (this 
corresponds to < a; < 1 in the previous notation). Thus the filled in Julia set is just 
the segment of the real axis with — 2 < 2; < 2. This analysis shows that the forward 
iterations send \z\ to infinity if the temperature has an imaginary component or if 
the magnetic field has both real and imaginary parts non-zero (this latter possibility 
would result in m having non-zero imaginary part and thus so would z). Thus chaotic 
trajectories occur only for real K and pure imaginary h. The behaviour of the Julia 
set under iteration is shown in Figure 2 for c = —2. 

An obvious question is: how generic is this behaviour? For a general Hamiltonian 
when is it possible to obtain chaotic behaviour in some region of (complex) coupling 
space? For the moment this question must remain unanswered, but a few comments 
should be made. Feigenbaum was aware of the universality in chaos [13] . Near an 
extremum any non-linear map (with non- vanishing second derivative) can be put into 
the form (13) with the number 4 replaced, in general, by a parameter, A. 

Thus 

x' = Xx{l - x) (22) 
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is generic, but whether or not one has chaotic behaviour depends on the value of A and 
the initial value of x. 1 do not know of any reason why A has the rather special value 
of 4 for the one dimensional Ising model. More generally, Feigenbaum has shown [13] 
that the properties of a general non-linear map 

x' = Xfix) (23) 

are independent of the exact form of f{x) near a maximum. For any particular model, 
the value of A would have to be calculated ab initio and I know of no way of deciding in 
advance whether or not chaotic flow would result. For the 1-D Ising model it is clear 
that the onset of chaotic trajectories is related to the second order phase transition 
at T = and the existence of the Lee- Yang edge singularity. However other models 
show chaotic recursive maps for real values of the couplings which are not related to 
second order phase transitions - rather they are related to frustration and glass-like 
structures [3]. Thus it does not appear that a second order phase transition is a 
pre-requisite for chaos, but neither are spin-glasses a pre-requisite. Unfortunately the 
number of models for which the recursion relations are known exactly is rather few. 
In particular the recursion relations for the 2-D Ising model are not known, so it is 
not possible at this stage to say whether or not the two dimensional Lee- Yang edge 
singularity is related to chaotic trajectories. 

It is a pleasure to thank Chris Stephens and Bill McGlinn for useful discussions on 
the one dimensional Ising model, Tony O'Farrell for helpful discussions about the 
Julia set and Paul Upton and Denjoe O'Connor for invaluable suggestions concerning 
Lee- Yang zeros. 
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Figure Captions 



Figure 1. Critical line in the K-9 plane. The renormalisation flow is regular below the critical 
line X — Q and chaotic above it for all values of •0 which are irrational (x = siv? {ttiI})) . 

Figure 2. A representation of the Julia set in the complex 2;-plane for c = —2. The Julia set 
itself is the real line segment —2 < 'R.e{z) < 2 which is the width of the diagram. The 
contours depict the rate at which a point is repelled from the Julia set, the darker 
the contour the less rapid the expulsion. The picture was generated use the program 
FRACTINT produced by the Stone Soup Group. 
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